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History 



► G.T. Fechner, et al. (1828-1900) 

► A. Lotka (1910-1920). 

► V. Volterra (1926) 

► B. Belousov (1951) 

► A.M. Zhabotinsky (1961) 

► Ilya Prigogine and his Brussels school 


BZ reaction-chemical reaction exhibiting oscillatory behavior. 
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Systems 

of Chemical Reaction 



► Open system 

Exchange of both molecules and energy with the surroundings 
is allowed, (for in vivo studies) 

► Closed system 

Exchange of energy but NOT molecules with the surroundings 
is allowed, (for in vitro studies) 


Most of currently exisiting reaction models exihiting oscillation are 
open system. 
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Irreversible Lotka-Volterra Model 


Irreversible Reaction. 


A + X^b 2X, X+Y^b2Y, Y B, 


► By the Law of Mass Action, one has 


Closed system 


Open system 


c A = — k\ c A x 


x = k\ c A x — k 2 xy 
y = k 2 xy - k 3 y 
, Cb = k 3 y. 


x = k\c A x — k 2 xy 
y = k 2 xy - k 3 y. 


( 1 ) 


(2) 


(3) 
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Reversible Lotka-Volterra System 


► Reversible Reaction. 

k\ l<2 

A + X^±2X, X+Y^±2Y, 

k - 1 k -2 


Y S B. (4) 

k-3 
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Reversible Lotka-Volterra System 


► Reversible Reaction. 


k\ 

A + X ^ 2X 

k - 1 

► Rate Equations 


dx 



dt 


l<2 /C 3 

X + Y — 2Y, Y 4 B. 

k- 2 k . 3 

kiCAX ~ /c_ix 2 - k 2 xy + /c_ 2 y 2 , 
k 2 xy - k- 2 y 2 - k 3 y + k~ 3 c B , 
-kiC A x + /c_ix 2 , 
k 3 y — k- 3 c B ■ 
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Nondimensionalization 



► Rescaling 

u = v = %y, w = %c A , z = |c B , t = fet, 

► Dimensionless Form. 


du 

= u(w — v) — e(au 2 — 

— 

dr 

dv 

= v(u — 1) — ev 2 + ez 

d T 

dw 

= —cr(wu — ecru 2 ). 

dr 

dz 


dr 

II 

1 

CO 

N 


(6) 


(7) 
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Closed System 



► Linear Conservation Law. 


w 


u + v -\ b z = £ = constant. 

<7 


Reduced system. 


du 


& 




dr 


u(w — v) — e(ou 2 — v 2 ) 

v{u - 1) - ev 2 + e (£ - u - v - 

—o(wu — sou 2 ). 


(8) 
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Closed System 



Denote 

T = j(tv, v. w) € M 3 , u,v,w > 0, and u + v + — < . 

Then T is positively invariant under the flow induced by the closed 
system (8), and T is called the reaction zone. 

System (8) has a unique interior equilibrium point P e T at which 
its Jacobian matrix has three real eigenvalues 

|Ai + (1 + e)| ~ e 2 , |A 2 + eC|~e 2 , |A 3 + cre 2 £\ ~ cr 2 e 3 . 

Thus P is an asymptotically stable node. 
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The Second Law of Thermodynamics 


P is the global attractor of system (8) in T. The free energy 


I - )+ vln(-) + -In + Inf ? " " “ * ' 


£ — u* — v* — — — 

’ CT 


serves as the Lyapunov function, where P = ( u *, v*, w*). 



< □ ► « ► < 1 ► < 1 ► = -OQ ,0 


Yongfeng Li (joint work with Yingfei Yi and Hong Qian) 


Nonlinear Oscillation and Multiscale Dynamics in a Closed Chi 




Outline 

Oscillatory Chemical Reaction 

Main Results 

Canonical vs Grand Canonical Systems 
Future Work 


Closed System - The Second Law of Thermodynamics 
Far-from-equilibrium Dynamics in T° 

Near-equilibrium Dynamics in T n 
Dynamical Transition from T 0 to T n 


Far-from-equilibrium - Oscillation Zone 


There exist T° C T and a ID curve \N° e C T° such that all the 
oscillatory solutions in T° go around W° e . The total oscillation 
time is T a — — — _ For fixed wq > w\, there are at least N 

complete oscillations by 2ir with w decreasing from wq to wi, 
where N > ln Wl — 1 . And W° e is approximately given by 
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Consider the partially perturbed system where e = 0 


u = u(w — v) 
v = v (u — 1 ) => 

w = 0 


which admits a stable invariant manifold 


w a = |(o. v, w), U = fJ, rT = - 

with Lyapunov function 

Sr = (! + *) 


cr 


u(w — v ) 

v(u-l) 

—awu 


(9) 


\/ = 1 / 1 / > 0 > . 


( u\ 


. / v M 

u- n a - In — 

+ 

1 

5 

1 

5 

3 

1 

V/W. 


V W / J 
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Generalized Lyapunov Type Numbers 


1l{W°) = 

a L (W°) = 


\^ dmw:w < 

\og\\ DMW2)^\\ 
-oo \og \\n^ D(j) t {W°)\\ 


lim 

t — y — oo 

lim 


1, 

> 2 , 


► Exponential Dichotomy 


u — > /J, a + x, v — > v + w, w — > w, 
(aX> = A„MX + G„ A X, W ) 
i W = F^^iX, w ), 

where Re(X(A a )) < — as w > 


X = (x,y) T 


0 

w 

f-^cr 


f^<7 

-Gila 
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Proof-Sakamoto (1990) 

1 



► Modified system 


f aX’ = A a (w)X + G a , £ {X,w) 

\ W = F CTj£ (X, w)x[ a 2^ ] {w) 

^ | x (t) = ^ j ^a{t,s,w{s))G a ^{X,w)ds 

\ w(t ) = oo, M/(0) = rj € [cr 2 ,^]. 

► .F(X) = 1 f <t> t7 (t,s,H(n,X)(s))G tT , e (X,H(ri,X)(s))ds. 

^ J —oo 


► J 7 is a contraction as |X| < 5 for some 5 => F{X*) = X*. 

► w° £ = {(u = // CT +x* ( 0 ), 1/ = w + y*( 0 ), 1/1/), W € [cr 2 ,<]}. 
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Closer Look of Oscillatory Behavior 


Define 

T° = {(tv, v,w) eT : w > a} , T 2 = {(tv, v, w) € T : w > a 2 } , 

= {(tv, v,w)£T : v < w} , Q + = {(tv, v, w) € T : v > w} . 

For e <C a <C 1 and some 1 < a < 2, 

► T° C T° C T 2 °. 

► In Q+, W2k+1 ~ 1 and ~ 1. 

► | n Q_ w 2fc+2 < e -(TCl(w 0 ,Eo) anc | £ 2*+2 < e -o-C 2 (wo,E 0 ) 

’ w 2k+l E 2k+1 

► At the bottom of T° where w ~ a a , E ~ (cr — 2)cr“ In a. 
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Near-equilibrium - Non-oscillation Zone 


There exist T n C T and a 2D strongly stable invariant manifold 
M" e C T n and a ID stable curve W" E C M" e . The “total time” 

in T n is T n ~ And 

w a,e ~ | (u, V, w), v = e w = aeu, u £ [0, u], u < 1 1 . 


Non-oscillation Zone 



x 10" 8 



Yongfeng Li (joint work with Yingfei Yi and Hong Qian) 


Nonlinear Oscillation and Multiscale Dynamics in a Closed Chi 




Outline 

Oscillatory Chemical Reaction 

Main Results 

Canonical vs Grand Canonical Systems 
Future Work 

Existence of M" 


Under the following transformation 

u — } u, v — >• ev , w — y crew. 

system (8) becomes 

= e [ u(aw — v) — (a u 2 — e 2 v 2 )] 

= v(u - 1) — e 2 v 2 + (£ — u - EV - ew) (12) 
= —au(w — u). 

By Fenichel's Theorem, critical manifold 

MS = j(u, v, w), v = u,w€ [0, D], u < 1 

is normally hyperbolic and thus persists under the perturbation. 


du 

fv 


d, 

d i 


dr 
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Existence of W" 



► Reduced System 


du 

dr\ 

dw 

dr\ 


- \u{aw - h ae {u, i/i/)) - (au 2 
a 

— u(w — u). 


Critical manifold 


- 2l2 


h l,e ( U > W ))] 


(13) 


Wq = {(u, w), w — u € [u, u], 0 < u < u < 1} 

is normally hyperbolic and thus persists under the 
perturbation. 
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In the Vicinity of Equilibrium 

P 


► Stable Invariant manifold may be applied. 

► Further Rescaling 


U — >• (T£U, v — >• V 

£ 

W — >■ W, T2 = T\. 

a 


yields 



-T ~ = u(aw - hve) - (C T 2 £U 2 - -h a £ ) 
dj 2 O 

dw 2 

— — = —a U{W — (T£U). 

OT2 

(14) 
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° to r n 


The passage of entering the non-oscillation zone T n from the 
oscillation zone T° is around the portion of the central axis 
connecting W° e and W" e . This is exactly where the transition 
occurs. 
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A General Network of Chemical Reactions 


Consider a system of chemical reactions whose rate equations, by 
the law of mass action, are given by 


X' = V{X) = AR{X ) 


where A = (a,y) is the stoichiometric matrix and 
Ri(X) = rf {X ) - r/»(X) with 


r[(X) = k, Xj a \ 

3/1 <0 

Rewrite equation (15) as 

X[ = 

X' = 


5 (x) = k_i n 


x j 


a , ; >0 


Fi(X 1 ,X 2 ) 
F 2 (X 1 ,X 2 ). 


(15) 


(16) 
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Canonical System and Grand Canonical System 


Canonical System. 

By the linear conservation law MiXi + M 2 X 2 = £, equation 
(16) is reduced into 

= F 1 {X 1 ,M^{X~M 1 X 1 )). (17) 

Grand Canonical System. 

By treating X 2 = X 2 ° as a constant vector, equation (16) is 
reduced into 

X[ = F 1 (X 1] X°). 


(18) 
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Canonical System and Grand Canonical System 


► (Gibb’s Principle) For given X 2 and £, let X* and X* c be the 
equilibrium points of systems (17) and (18), respectively. If 
X 2 ° = (f - M x X*), then X* = X* c . 

► Systems (17) and (18) share the "same” Lyapuov functions. 

► Both X* and X* are all asymptotically stable nodes. 
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Recall 


du 

I 

t 

& 

dr 


Canonical 


du 

I 

qf T 

dw 
, dr 


u(w — v) — e(au 2 — v 2 ) 
v(u — 1 ) — ev 2 + ez 
—a(wu — ecru 2 ) 
v — ez. 


u(w — v) — e(<ju 2 — v 2 ) 
v( t l-l)- £ V 2 + e (*- U -V-?) 
—a(wu — sau 2 ) 


Grand Canonical 


du 

fv 



u(w — v) — e(<ju 2 — V 2 ) 
v(u — 1 ) — EV 2 + EZ 
—a(wu — Ecru 2 ) 

< n ► < S 1 ► < = 


( 19 ) 


( 20 ) 


( 21 ) 
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Canonical System vs Grand Canonical System 


System (20) admits similar dynamics as (21) does. 
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Canonical System vs Grand Canonical System 


► When v = v* ~ s. 
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Future Work 


► Effect of noise on the dynamics 

► Macroscopic level- Fokker-Planck equation 

► Microscopic level- chemical master equation(in progress) 

► Dissipative perturbation of conserved system. 
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